Primordial black hole formation from non-Gaussian 

curvature perturbations 



p. A. Klimai'^: E. V. Bugaev"^ 
Institute for Nuclear Research, Russian Academy of Sciences, 
60th October Anniversary Prospect la, 117312 Moscow, Russia 



Abstract 

We consider several early Universe models that allow for production of large curvature 
perturbations at small scales. As is well known, such perturbations can lead to formation of 
primordial black holes (PBHs). We briefly review the today's situation with PBH constraints 
and then focus on two models in which strongly non-Gaussian curvature perturbations are 
predicted: the hybrid inflation waterfall model and the curvaton model. We show that PBH 
constraints on the values of curvature perturbation power spectrum amplitude are strongly 
dependent on the shape of perturbations and can significantly (by two orders of magnitude) 
deviate from the usual Gaussian limit Vc^ < 10~^. We give examples of PBH mass spectra 
calculations for both inflationary models. 

U; 

p I 1 Introduction 

As is well known, in models of slow-roll inflation with one scalar field the curvature per- 
turbation originates from the vacuum fluctuations during inflationary expansion, and these 
fluctuations lead to practically Gaussian classical curvature perturbations with an almost flat 
power spectrum. However, it is well known also that both these features are not generic in the 
case of inflationary models with two (or more) scalar fields: such models can easily predict adia- 
batic perturbations with, e.g., a "blue" spectrum and these perturbations can be non-Gaussian 

Possibilities for appearing of non-Gaussian fluctuations in inflationary models with multiple 
scalar fields had been discussed long ago [21 El H]. The time evolution of the curvature pertur- 
bation on super horizon scales (which is allowed in multiple-field scenarios [5]) implies that, in 
principle, a rather large non-Gaussian signal can be generated during inflation. According to 
the observational data [6], the primordial curvature perturbation is Gaussian with an almost 
scale-independent power spectrum. So far there is only a weak indication of possible primordial 
non-Gaussianity [at (2 — 3)(T level] from the cosmic microwave background (CMB) temperature 
information data (see, e.g., [7j). However, non-Gaussianity is expected to become an important 
probe of both the early and the late Universe in the coming years [8]. 

The second important feature of predictions of two-field models is that these models can 
lead to primordial curvature perturbations with blue spectrum (for scales which are smaller 
than cosmological ones) and, correspondingly, can predict the primordial black hole (PBH) 
production at some time after inflation. In this case, PBHs become a probe for the non- 
Gaussianity of cosmological perturbations [9l [lOl [TTl [12]. The results of PBH searches can be 
used to constrain the ranges of early Universe model parameters. 

There are several types of two-field inflation scenarios in which detectable non-Gaussianity 
of the curvature perturbation can be generated: curvaton models [131 E US US ES], models 
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with a non-inflaton field causing inhomogeneous reheating |17^ I18j . curvaton-type models of 
preheating (see, e.g., [19] and references therein), models of waterfall transition that ends the 
hybrid inflation [20l [211 [211 [23] . In these two-field models the primordial curvature perturbation 
^ has two components: a contribution of the inflaton (almost Gaussian) and a contribution of 
the extra field. This second component is parameterized by the following way [24j 

G(x)=aa(x) + a2(x)-(c72). (1) 

If a = 0, one has a x^-model. Obviously, the quadratic term can't dominate in ^ on cosmological 
scales where CMB data are available. It can, however, be important on smaller scales. 

In the present work we study the predictions of the PBH production for two particular two- 
field models: the hybrid inflation model with tachyonic instability at the end of inflation and the 
curvaton model. The potentially large non-Gaussianity in these models is connected with the 
fact that the predicted magnitude of the curvature perturbation is proportional to a square of 
the non-inflaton (waterfall or curvaton) field. The blue spectrum in the hybrid waterfall model 
can arise, e.g., through the tachyonic amplification due to the dynamical symmetry breaking 
[20] and, in the curvaton model, due to, e.g., supergravity effects leading to the large effective 
mass of the curvaton [Ij. 

The main attention in the present paper is paid to a study of probability distribution 
function (PDF) of the curvature perturbation and the shape of the black hole mass function, 
with taking into account of the non-Gaussianity. The first general study of PDF of the curvature 
perturbation in curvaton model was carried out in |25j . 

Primordial curvature perturbation spectrum in hybrid inflation waterfall model was recently 
studied in[26l[271[28l[29l[30l[3ll[32l[33l[31], PBH production in this model was considered in [32l 
[33l [Ml [35] . PBH production in curvaton scenario was studied in recent works [36l [37] (without 
considering the non-Gaussian effects). The approximate PBH constraints on the curvature 
perturbation power spectrum in the curvaton model were obtained in |i34j . 

The plan of the paper is as follows. In Sec. [21 we briefly review the available constraints 
on PBH abundance that follow from different types of astrophysical and cosmological data. In 
Sec. [3|we present the recent limits on the amplitudes of curvature perturbation power spectrum 
obtained with an assumption of Gaussianity of primordial curvature perturbations. In Sec. [Dthe 
hybrid inflation model and constraints on its parameters, coming from PBH non-observation, 
are discussed. In Sec. [5] we discuss the possible production of PBHs in the curvaton model and 
the corresponding cosmological constraints that can be obtained. Our conclusions are given in 
Sec. M 

2 Available constraints on PBH abundance 

There are several sources of information that allow to obtain limits on PBH abundance. In 
the region of Mbh which is of most interest for us (10^ < Mbh ^ 10^^ g) these limits can be 
divided in three groups: 

(i) constraints on PBHs from big bang nucleosynthesis (due to hadron injections by PBHs 
[38], photodissociation of deuterium [39] and light nuclei, fragmentations of quarks and gluons 
evaporated by PBHs t40j), 10^ < Mbh ^ 10^^ g, and from influence of PBH evaporations on 
the CMB anisotropy, 2.5 x lO^^ < Mbh < 2.5 x 10^'' g [1I|, 

(ii) constraints on PBHs from extragalactic photon background [12], 10^'^ < Mbh ^ 10^'' g, 

(iii) constraints on non-evaporating PBHs (gravitational and lensing constraints, Mbh > 
10^5 g). 

Constraints on PBHs from data on extragalactic neutrino background |431l44j. in the region 
10^^ < Mbh ^ 10^^ g, are somewhat weaker than nucleosynthesis constraints. In the PBH mass 
region Mbh ~ 10^^ g there also exist constraints from non-observation of induced gravitational 
waves, which have been studied in [451 146|. 
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Figure 1: The available limits on PBH abundance /3'{Mbh) from different types of observations. 
The figure is taken from 



For a detailed review of these and other PBH constraints, see |4H I47j. Traditionally, the 
constraints are expressed as limits on energy density fraction of the Universe contained in 
PBHs at the moment of their formation, Ppbh^ as a function of PBH mass Mbh (it is typically 
assumed that 

Mbh = fhMh, (2) 

where is the horizon mass at the moment of PBH formation and fh is a constant of order 
of one). The summary of existing PBH limits, taken from [5T], is shown in Fig. [TJ This Figure 
gives limits in terms of /3', which is simply related to /3pbh j41j: 

P'{Mbh) = fl'' {-^y' Ppbh{Mbh), (3) 

where g.^ is the effective number of relativistic degrees of freedom. 

The available limits on PBH abundance, shown in the figure, can be transformed in con- 
straints on other cosmological parameters. If, for example, inflationary model predicts large 
values for the curvature perturbation power spectrum at some cosmological scale, the model 
parameters can be constrained based on the requirement that PBHs are not over-produced. 



3 PBH constraints in the Gaussian case 

Several single-field inflationary models which predict a peaked power spectrum at some 
scale ^0 have been considered in [48 [ I49j . The corresponding constraints on the amplitudes 
of the curvature perturbation power spectrum from the non-observation of PBHs have been 
obtained in [36], with an assumption of Gaussian form of the perturbation distributions. It 
had been shown in [3^ for the particular case of Coleman- Weinberg inflationary potential that 
non-Gaussianity of perturbations is, really, rather small. Using the assumption of Gaussianity, 
the PBH abundances can be calculated with a probability distribution of the form 



2al(M) 
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Figure 2: PBH bounds on the maximal value of primordial power spectrum, "P^, for the 
Gaussian case. Constraints depend on the peak's width S and its position (which corresponds 
to PBH mass Mbh)- The figure is taken from |46j . 



where aji{M) is the mass variance (mean square deviation of density contrast 6 in the sphere 
of comoving radius R), M is the initial mass of fluctuation. 

It is convenient to use some kind of parametrization to model the realistic peaked power 
spectrum of finite width. In [l6] the distribution of the form 

lgVn{k) = B + {lgVl-B)e^p 

was used. Here, B ~ —8.6, characterizes the height of the peak, is the position of the 
maximum and S is the peak's width. Note that comoving curvature perturbation TZ practically 
coincides with uniform density curvature perturbation on super-horizon scales. 

Calculating the abundances of PBHs for different sets of model parameters, and using the 
available PBH limits [H], the constraints on the value of for different values of S were 
obtained in [IB]. The corresponding results are shown in Fig. [5J It is seen from this Figure that, 
with the assumption of Gaussianity, the resulting constraints on the power spectrum amplitude 
are of the order of 10~^ with some scale dependence. Similar results were also obtained recently 
in 07]. 




4 The waterfall transition model 

In this section we consider the hybrid inflation model which describes an evolution of the 
slowly rolling inflaton field (f) and the waterfall field X: with the potential [50^ [5T] 

Vi<P, x) = [m^ - ^X"^ + \m'<t^' + \ict>\'. (6) 

The first term in Eq. ([G]) is a potential for the waterfall field x with the false vacuum at x = 
and true vacuum at Xq = 2M^/\/A = v'^ . The effective mass of the waterfall field in the false 
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Figure 3: The numerically calculated spectra V^^k) at the moment of the end of the waterfall. 
Parameters used for the calculation are: r = 0.1, He = 10^^ GeV (all curves); </>c = 3 x W~^Mp 
(for /? = 1, 2, 7); (/)c = O.lMp (for /3 = 100, 1500). The figure is taken from [33j. 



vacuum state is given by 



m2(</,) = 7(</,2_02)^ 
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(7) 



At 0^ > the false vacuum is stable, while at 0^ < (j)'^ the effective mass-squared of x becomes 
negative, and there is a tachyonic instability leading to a rapid growth of x-niodes and eventually 
to an end of the inflationary expansion [30j • It is convenient to define the following parameters 



c 




(8) 



where He is the Hubble rate during inflation. 

It was shown in [33J that the magnitude of the curvature perturbation can be presented in 
the form 

C = Cx = -^(x'-(x')), (9) 

where and (x^) are determined at the time of an end of the waterfall, t = tend, and A is 
given by the integral 



Hedt 



fit) 



(10) 



Here, the function f{t) describes the time evolution of the waterfall field (it is proportional 
to a numerically obtained solution of a linearized field equation for a Fourier component of a 
waterfall field perturbation, Sxk{t), which is almost independent on A;). Energy density p and 
pressure p are a sum of the contributions of (j) and x fields. 

It was shown in [33] that for /3 ~ 1 the curvature perturbation spectrum will reach values of 
"P^ ~ 1 in a broad interval of other model parameters (such as r, 7, He). The peak values, /c*, 
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for small /3, are far beyond horizon, so, the smoothing over the horizon size will not decrease 
the peak values of the smoothed spectrum. The examples of the calculated power spectra for 
this model are shown in Fig. [3l 

The relation between curvature perturbation and the waterfall field value is given by Eq. 
dSD, or, using ^ = (x^), 

( = -A{x^-al) = Cmax-Ax^ Cmax = Aal. (11) 

Here, A and o"^ generally depend on the smoothing scale R. The distribution of x is assumed 
to be Gaussian, i.e., 

1 -4 

P^(X) = 7=e ^-x. (12) 

o"-^ V zvr 

The distribution of C can be easily obtained from (|lll I12p : 

dx 



PdO =Px 



dC 



\/'^'^Cmax {Cmax C) 



e 2Cma:r , ( < (^^^ ^ (13) 



which is just a ^^-distribution with one degree of freedom, with an opposite sign of the argument, 
shifted to a value of (max- As required, {0 = and 



Cmax 

>2 



On the other hand, one has 



{C)= / CpdOdC = 2C^,,. (14) 



(0=^C= I V^{k)W\kR)^, (15) 



where W{kR) is the Fourier transform of the window function, and we use a Gaussian one, 
W'^{kR) = exp(— /c^i?^), in this work. 

From (fT^ we can write for (max (we now denote the argument R explicitly) : 



Cmax (-^) 



i J V,{k)W\kR)'^ 



ni/2 

(16) 



Below, we will use the following notation: Cmax{R = 0) = Cmax- So, 



1 n.dk' 

2 J P^^^-'T 



1/2 , 

= ^icy^'- (17) 



It is clear that PBHs can be produced in the early Universe, if Cmax > Cc where Cc is the thresh- 
old of PBH formation in the radiation-dominated epoch, which is to be taken from gravitational 
collapse model. For estimates, in the following we will use two values: Cc = 0.75 and Cc = 1 
(corresponding to the PBH formation criterion in the radiation-dominated epoch: 5 > 6c, with 
(5c = 1/3 and 6c = 0.45, respectively) [35] . 

The energy density fraction of the Universe contained in collapsed objects of initial mass 
larger than M in Press-Schechter formalism [52] is given by 



oo 



- I Mn{M)dM = j p^{C)dC = P{C > Cc;RiM),U), (18) 

where function P in right-hand side is the probability that in the region of comoving size R 
the smoothed value of C will be larger than the PBH formation threshold value, n{M) is the 
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mass spectrum of the collapsed objects, and pi is the initial energy density. Here we ignore the 
dependence of the curvature perturbation C, on time after the end of the waterfall, assuming it 
does not change in super-horizon regime, until the perturbations enter horizon aX, k = aH. 

The horizon mass corresponding to the time when fluctuation with initial mass M crosses 
horizon is (see [44j ) 

Mh = m//^M2/3, (19) 
where Mi is the horizon mass at the moment ti, 

M, « —t\p, p. —{H~^fp = (20) 

(here, we used Friedmann equation, pi = 3MpH^). The reheating temperature of the Universe 
is 



( 



90Mf,H, 



1/4 



Trh = 2^ , 5* ~ 100. 21 

For simplicity, we will use the approximation that mass of the produced black hole is pro- 
portional to horizon mass [see Eq. namely, 

Mbh = fhMh = hMl'^'M^I^ (22) 

where fh ~ (l/S)-*^/^ = const. 

Using (jlSp and ()22p . for the PBH number density (mass spectrum) one obtains |35j 



'47r\-^/=^ 
UBHiMBu) = ( Y ) 



dP 



OR 



(23) 



Considering the moment of time for which horizon mass is equal to M/j, one can obtain for 
the energy density fraction of Universe contained in PBHs: 

1 / Mh \ ^^"^ f 

^PBH{Mh)- j UBHMlHdlnMBH- (24) 

It is well known that for an almost monochromatic PBH mass spectrum, QpBHiMh) coincides 
with the traditionally used parameter Ppbh- Although all PBHs do not form at the same 
moment of time, it is convenient to use the combination p^ M^^nBHiMBn) to have a 



feeling of how many PBHs actually form, i.e., to use the estimate following from 

Mr'/^pr^M'J^nBH{MBH) « ^pbh- (25) 

Fig. m shows the examples of PBH mass distributions produced in the waterfall model, 
calculated assuming the curvature perturbation power spectrum has the form 



(Igk/kof 



2S2 



(26) 



where gives the maximum value approached by the spectrum, ko is the comoving wave 
number corresponding to the position of the maximum and S determines the width of the 
spectrum (this is just a parametrization of the result of Fig. [3D. 

It is seen from Fig. S] that in the waterfall model PBH abundance severely depends on the 
amplitude of the curvature perturbation spectrum in a fine-tuning regime: once (max is above 

PBHs are produced intensively. Demanding that no excessive amount of PBHs forms in the 
early Universe, we can impose the bound on parameters of the inflaton potential. From the 
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Figure 4: The PBH mass spectra for different values of the perturbation spectrum ampHtudes 
for the hybrid inflation waterfall model. From right to left, = l{Cmax ~ 1-42), 0.4(^maa; ~ 
0.9),0.28(Cmax ~ 0.752), 0.27846(Cmaa; ~ 0.750012), 0.2784511(Cmax ~ 0.750000068). The posi- 
tion of the peak in 'P(^(A:)-spectrum is the same for all cases. For the calculation we used the 
value E = 0.7, and Cc = 0.75. The mass Mj^ corresponds to horizon mass at the moment of 
time when perturbation with comoving wave number /cq enters horizon. The figure is taken 
from [35]. 

condition Cmax < Cc one has, for two fixed values of Co the following constraints (r = 0.1; there 
is a weak dependence on this parameter but the result is almost independent on 7, He): 

Cc = 0.75 : /3 > 2.3, < 0.29; (27) 

Cc = 1 : /3 > 1.65, < 0.55. (28) 

One can see that the limit on the value of V,^ is significantly weaker than in the Gaussian case 
(see Fig. [2]). 

5 A case of the curvaton-type model 

Curvaton is an additional to the infiaton scalar field that can be responsible (partly or fully) 
for generation of primordial curvature perturbations [U [T^ \T5\ [16] . This field can also be a 
source of PBHs, as discussed in [53l [32]. 

In this work, we only consider the case of a strong positive tilt (the possibility discussed 
in [H [16]) of the curvaton-generated perturbation power spectrum. At the same time, it is 
assumed that infiaton is responsible for generation of perturbations on cosmological scales (see 
Fig. [5] for an illustration). 

The curvaton field generates cosmological perturbations in two stages [Tl [H] [T5 ] [16]: 

(i) Quantum fiuctuations of the curvaton during infiation (at time of horizon exit) become 
classical, super-horizon perturbations. 

(ii) In the radiation-dominated stage, the curvaton starts to oscillate (this happens at the 
time when Hubble parameter becomes of order of curvaton's effective mass, H ~ m). The 
Universe at this stage becomes a mixture of radiation and matter (the curvaton behaves as 
a non-relativistic matter in this regime). The pressure perturbation of this mixture is non- 
adiabatic and the curvature perturbation is thus generated. One obtains, approximately, the 
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Figure 5: A sketch that ihustrates a relation between curvaton-generated and inflaton- 
generated curvature perturbation power spectra for the scenario of PBH production that we 
consider. 



expression (see, e.g., [53] ) 



C(i,x) 



rcr' r / fj' 



2(7n 



(29) 



where r is the density parameter, r = Apfj/{Apr + 3/9^) (/Or is the energy density of radiation 
after inflation), dose is the value of the curvaton field at the onset of oscillations. The initial 
value for the curvaton field, 5cj*, is set by inflation. The derivative in Eq. (|29p is taken with 
respect to the field value during inflation, cr*. The term containing the second derivative, a'^^^^ 
is neglected. It is assumed that r <C 1. 

Assuming zero average value for the curvaton field (i.e., working with the maximal box |16j). 
we keep in Eq. ([29j) only the second term. 



C(t,x) = ^f^V5a2 

4 V CJn, 



(30) 



The fluctuations are strongly non-Gaussian which is not forbidden on small scales. Note that 
the sign of the perturbation in Eq. (|30|) is different from the one in Eq. @. 

The curvaton-generated curvature perturbation spectrum can be written [16] as (using the 
Bunch-Davies probability distribution for the perturbations of the curvaton field [551 [56]) 



V. 



1/2 



Ml" f 



k 



(31) 



where VL^j = / p is the relative curvaton energy density at the time of its decay (it depends on 
the curvaton decay rate) and 



Here, is the curvaton effective mass during inflation and H^, is the Hubble parameter. 

It is rather natural (see [1] and, e.g., [57], which considers the models of chaotic inflation 
in supergravity) that ~ 2/3 which corresponds to a blue perturbation spectrum with the 
spectral index 

n = l + 2t^r^ 7/3 (33) 
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Figure 6: The form of the distribution (j36|) for cr^ = 2x 10 ^. Dashed lines show the considered 
values of Cc (0.75 and 1, see Sec. S]). 



(such a situation is shown in Fig. [5]). For the following, we parameterize the spectrum (|3ip in 
a simple form 



n-l 



(34) 



and will treat V^, n, and fcg as free parameters. Note that does not, in general, coincide 
with the comoving wave number corresponding to the end of inflation. Rather, it is the scale 
entering horizon at the time when is created |16j . 

Using (j30p and subtracting the average value, we can write for the curvature perturbation 
the formula 

C = A{6a'-{6a')), (35) 
and, in analogy with the case of Eq. (|13|) . for the distribution of C one obtains 



pdO 



C>Cn 



where 



Note also that 



-A{5a^) < 0. 



(36) 



(37) 



(38) 



The form of the distribution 



for cr^ 



2 X 10 ^ is shown in Fig. It is seen that in this 



particular case the probability to reach ~ ~ 1 is 



10 



-20 
-4 ;c 



i.e., of the same order as PBH 



constraints on /3' of Fig. m Thus, the value of V^{k) ~ lO""^ is aheady enough to produce an 
observable amount of PBHs in this model (this is in agreement with the estimates of j32t 134)). 

For a more exact estimate, we can calculate PBH mass distributions that are generated for 
a particular set of parameters (n, V^, etc.) and then compare the resulting f^pBH (using Eq. 
(p5|l ) with the known limits [IT] (see Fig. [T|). 

The example of PBH mass spectrum calculation is shown in Fig. [71 The formula ()23p from 
the previous Section is used for this calculation (using of Eq. ([36]) for the calculation of the 
probability P). It is seen from Fig. [7] that PBH abundances depend on particular choice of Cc- 

The resulting constraints on parameter (for two values of the spectral index n) from 
PBH non-observation, that were obtained in this work, are shown in Fig. [51 Thus, we see 
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Figure 7: Examples of PBH mass spectra calculated for the curvaton model. The parameters 
used are: n = 2, = A x 10"'*. For the upper curve, Cc = 0.75, for the lower one, Cc = 1- The 
mass Mj^ corresponds to horizon mass at the moment of time when perturbation with comoving 
wave number kg enters horizon. 

that the available constraints on the PBH abundance can be transformed into limits on the 
power spectrum amplitude and, subsequently, on curvaton model's parameters, such as Q^- For 
example, comparing Fig. [8] and Eq. ()3ip . one obtains, roughly, 

« (^C)'^' ^ 10"'' (39) 

which is already a useful constraint. It is a subject of our further study to get more exact limits 
for particular parameter sets of the model. 

logjo Mh [g] 
30. 28. 26. 24. 22. 20. 18. 16. 14. 12. 




10 12 14 16 
logjo ke [Mpc"^] 

Figure 8: The limits on the maximum value of curvature perturbation power spectrum 
from PBH non-observation, for the curvaton model. The forbidden regions are above of the 
curves. is the horizon mass at the moment of time when the perturbation with comoving 
wave number enters horizon. The spectrum is parameterized as Cc=0.75. 
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6 Conclusions 



Primordial black holes can be used to probe perturbations in our Universe at very small 
scales, as well as other problems of physics of early stages of the cosmological evolution. We 
have considered the PBH formation from primordial curvature perturbations produced in two 
cosmological models: in the waterfall transition in hybrid inflation scenario and in the curvaton 
model. Both models which were considered in the paper predict the production of strongly 
non-Gaussian perturbations, and non-Gaussianity was taken into account in the calculation of 
PBH mass spectrum (in Press-Schechter formalism). 

For both models, limits on the values of perturbation power spectrum as well as constraints 
on inflation model parameters were obtained. The approximate constraints on the curvature 
perturbation spectrum amplitude follow from Eqs. ()27l [28|) in the case of hybrid inflation 
waterfall model and from Fig. [8] in the case of the curvaton models. One can conclude that 
the constraint on "P^ in the case of hybrid inflation waterfall model is very weak, 'P(^ ^ 1. 
The corresponding constraint in the curvaton model is much stronger, "P^ < (10~^ — 10^'^), 
depending on the value of the spectral index and PBH mass. The PBH constraints obtained in 
this work confirm the estimates given in \30\ \3i\ . 



Acknowledgments 

The work was supported by The Ministry of education and science of Russia, project No. 
8525. 

References 



[1 

[2' 

[4 

[5] 
[6] 



[8] 
[9] 

[lo: 
[11 

[12 
[13 

[14 



A. D. Linde and V. F. Mukhanov, Phys. Rev. D 56 (1997) 535 [astro-ph/9610219] . 
D. S. Salopek, J. R. Bond and J. M. Bardeen, Phys. Rev. D 40 (1989) 1753. 

D. S. Salopek, Phys. Rev. D 45 (1992) 1139. 

Z. H. Fan and J. M. Bardeen, "Predictions of a nonGaussian model for large scale struc- 
ture," preprint UW-PT-92-11 (1992). 

A. A. Starobinsky, JETP Lett. 42 (1985) 152 [Pisma Zh. Eksp. Teor. Fiz. 42 (1985) 124 ]. 

WMAP Collab. (E. Komatsu et al), Astrophys. J. Suppl. 192 (2011) 18 [a rXiv:1001. 4538] 
[astro-ph.CO]]. 

A. P. S. Yadav and B. D. Wandelt, Phys. Rev. Lett. 100 (2008) 181301 [a rXiv:0712.lT48] 
[astro-ph]]. 

E. Komatsu, N. Afshordi, N. Bartolo et al. larXiv:0902.4759l [astro-ph.CO]. 

J. S. Bullock and J. R. Primack, Phys. Rev. D 55 (1997) 7423 |astro-ph/9611106] . 



P. Ivanov, Phys. Rev. D 57 (1998) 7145 [astro-ph/9708224| 



P. Pina Avelino, Phys. Rev. D 72 (2005) 124004 |astro-ph/0510052] . 
J. C. Hidalgo, larXiv:0708.3875 ^ [astro-ph]. 
S. Mollerach, Phys. Rev. D 42 (1990) 313. 

D. H. Lyth and D. Wands, Phys. Lett. B 524 (2002) 5 [hep-ph/0110002] . 



12 



[15] T. Moroi and T. Takahashi, Phys. Lett. B 522 (2001) 215 [Erratum-ibid. B 539 (2002) 
303] (hep-ph/0110096|. 

D. H. Lyth, J. Cosmol. Astropart. Phys. 0606 (2006) 015 | iastro-ph/06022 85|. 

G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69 (2004) 023505 
[arXiv: astro-ph /030359 1 1 . 



L. Kofman, arXiv:astro-ph/0303614 



K. Kohri, D. H. Lyth and C. A. Valenzuela- Toledo, J. Cosmol. Astropart. Phys. 1002 
(2010) 023 [Erratum-ibid. 1009 (2011) EOl] [a rXiv:0904.0793l [hep-ph]]. 

G. N. Felder, J. Garcia-Bellido, P. B. Greene, L. Kofman, A. D. Linde and I. Tkachev, 
Phys. Rev. Lett. 87 (2001) 011601 [hep-ph/0012142j . 



T. Asaka, W. Buchmuller and L. Covi, Phys. Lett. B 510 (2001) 271 [hep-ph/ 0104037]. 

E. J. Copeland, S. Pascoli and A. Rajantie, Phys. Rev. D 65 (2002) 103517 
[hep-ph/0202031|. 

J. Garcia-Bellido, M. Garcia Perez and A. Gonzalez- Arroyo, Phys. Rev. D 67 (2003) 103501 
[hep-ph /0208228]. 

L. Boubekeur and D. H. Lyth, Phys. Rev. D 73 (2006) 021301 |astro- ph/0504046]. 



M. Sasaki, J. Valiviita and D. Wands, Phys. Rev. D 74 (2006) 103003 [astro-ph/0607627] . 
D. H. Lyth, larXiv: 1005.24611 

J. Fonseca, M. Sasaki and D. Wands, J. Cosmol. Astropart. Phys. 1009 (2010) 012 
[arXiv:1005.4053j . 

A. A. Abolhasani and H. Firouzjahi, Phys. Rev. D 83 (2011) 063513 [arXiv: 1005.2934] . 
J. O. Gong and M. Sasaki, J. Cosmol. Astropart. Phys. 1103 (2011) 028 [arXiv: 1010.3405] . 



D. H. Lyth, J. Cosmol. Astropart. Phys. 1107 (2011) 035 [ arXiv:1012.4617 [astro-ph.CO]]. 

A. A. Abolhasani, H. Firouzjahi and M. Sasaki, J. Cosmol. Astropart. Phys. 1110 (2011) 
015 [ar Xiv:1106.6315 [astro-ph.CO]]. 

D. H. Lvth. larXiv:1107.1681l [astro-ph.CO], 

E. Bugaev and P. Klimai, J. Cosmol. Astropart. Phys. 1111 (2011) 028 [a rXiv:1107.3754] 
[astro-ph.CO]]. 

D. H. Lyth, J. Cosmol. Astropart. Phys. 1205 (2012) 022 [ arXiv:1201.43T? [astro-ph.CO]]. 

E. Bugaev and P. Klimai, Phys. Rev. D 85 (2012) 103504 |arXiv:1112.,560T [astro-ph.CO]]. 

M. Kawasaki, N. Kitajima and T. T. Yanagida. [arXiv: 1207.2550] [hep-ph]. 

H. Firouzjahi, A. Green, K. Malik and M. Zarei. [ajXiv:1209.2652l [astro-ph.CO]. 

Ya. B. Zeldovich, A. A. Starobinsky, M. Yu. Khlopov and V. M. Chechetkin, Sov. Astron. 
Lett. 3 (1977) 110. 

[39] D. Lindley, Mon. Not. R. Astron. Soc. 193 (1980) 593. 



13 



[40] K. Kohri and J. Yokoyama, Phys. Rev. D 61 (2000) 023501 far Xiv:astro-ph/9908160] . 

[41] B. J. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, Phys. Rev. D 81 (2010) 104019 
[arXiv:0912.5297 [astro-ph.CO]]. 

[42] D. N. Page and S. W. Hawking, Astrophys. J . 206 (1976) 1. 

[43] E. V. Bugaev and K. V. Konishchev, Phys. Rev. D 66 (2002) 084004 
[arXiv:astro-ph/02 06082| . 

[44] E. Bugaev and P. Klimai, Phys. Rev. D 79 (2009) 103511 jarXiv:0812.4247] [astro-ph]]. 

[45] R. Saito and J. Yokoyama, Phys. Rev. Lett. 102 (2009) 161101 |arXiv:0812.4339] [astro-ph]] . 



[46] E. Bugaev and P. Klimai, Phys. Rev. D 83 (2011) 083521 [ arXiv:1012.4697, [astro-ph.CO]]. 

[47] A. S. Josan, A. M. Green and K. A. Malik, Phys. Rev. D 79 (2009) 103520 [ arXiv: 0903. 3184) 
[astro-ph.CO]]. 

[48] R. Saito, J. Yokoyama and R. Nagata, J. Cosmol. Astropart. Phys. 0806 (2008) 024 
[arX iv:0804. 3470 [astro-ph]]. 

[49] E. Bugaev and P. Khmai, Phys. Rev. D 78 (2008) 063515 jarXiv:0806.4541 ' [astro-ph]]. 

[50] A. D. Linde, Phys. Lett. B 259 (1991) 38. 



[51] A. D. Linde, Phys. Rev. D 49 (1994) 748 [astro-p h/9307002| . 
[52] W. H. Press and P. Schechter, Astrophys. J. 187 (1974) 425. 

[53] K. Kohri, D. H. Lyth and A. Melchiorri, J. Cosmol. Astropart. Phys. 0804 (2008) 038 
[ar Xiv:0711.500 6 [hep-ph]]. 



[54] K. Enqvist and S. Nurmi, J. Cosmol. Astropart. Phys. 0510 (2005) 013 |; astro-ph/0508573] . 

[55] T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. Lond. A 360 (1978) 117; A. Vilenkin 
and L. H. Ford, Phys. Rev. D 26 (1982) 1231; A. D. Linde, Phys. Lett. B 116 (1982) 335. 

[56] A. A. Starobinsky, Phys. Lett. B 117 (1982) 175. 

[57] V. Demozzi, A. Linde and V. Mukhanov, J. Cosmol. Astropart. Phys. 1104 (2011) 013 
[arXiv: 1012.05491 [hep-th]]. 



14 



